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Abstract 

An elementary introduction is given to the subject of Supersymmetry in Quantum 
Mechanics which can be understood and appreciated by any one who has taken a first 
course in quantum mechanics. We demonstrate with explicit examples that given a solv- 
able problem in quantum mechanics with n bound states, one can construct new exactly 
solvable n Hamiltonians having n-l,n-2,...,0 bound states. The relationship between the 
eigenvalues, eigcnfunctions and scattering matrix of the supersymmetric partner poten- 
tials is derived and a class of reflectionless potentials are explicitly constructed. We 
extend the operator method of solving the one-dimensional harmonic oscillator prob- 
lem to a class of potentials called shape invariant potentials. It is worth emphasizing 
that this class includes almost all the solvable problems that are found in the standard 
text books on quantum mechanics. Further, we show that given any potential with 
at least one bound state, one can very easily construct one continuous parameter fam- 
ily of potentials having same eigenvalues and s-matrix. The supersymmetry inspired 
WKB approximation (SWKB) is also discussed and it is shown that unlike the usual 
WKB, the lowest order SWKB approximation is exact for the shape invariant poten- 
tials and further, this approximation is not only exact for large quantum numbers but 
by construction, it is also exact for the ground state. Finally, we also construct new 
exactly solvable periodic potentials by using the machinery of supersymmetric quantum 
mechanics. 



1 Introduction 



Supersymmetry (SUSY) is a symmetry between fermions and bosons. It was first intro- 
duced in High Energy Physics in an attempt to obtain a unified description of all basic 
interactions of nature. It is a highly unusual symmetry since fermions and bosons have 
very different properties. For example, while identical bosons condense, in view of Pauli 
exclusion principle, no two identical fermions can occupy the same state! Thus it is quite 
remarkable that one could implement such a symmetry. The algebra involved in SUSY is a 
graded Lie algebra which closes under a combination of commutation and anti-commutation 
relations. In the context of particle physics, SUSY predicts that corresponding to every ba- 
sic constituent of nature, there should be a supersymmetric partner with spin differing by 
half-integral unit. Further it predicts that the two supersymmetric partners must have iden- 
tical mass in case supersymmetry is a good symmetry of nature. In the context of unified 
theory for the basic interactions of nature, supersymmetry predicts the existence of SUSY 
partners of all the basic constituents of nature, i.e. SUSY partners of 6 quarks, 6 leptons 
and the corresponding gauge quanta (8 gluons, photon, Z°,W ). The fact that no scalar 
electron has been experimentally observed with mass less than about 100 GeV (while the 
electron mass is only 0.5 MeV) means that SUSY must be a badly broken symmetry of 
nature. Once this realization came, people started to understand the difficult question of 
spontaneous SUSY breaking in quantum field theories. It is in this context that WittenpQ 
suggested in 1981 that perhaps one should first understand the question of SUSY breaking 
in the simpler setting of nonrelativistic quantum mechanics and this is how the area of 
SUSY quantum mechanics was born. 

Once people started studying various aspects of supersymmetric quantum mechanics 
(SQM), it was soon clear that this field was interesting in its own right, not just as a model 
for testing concepts of SUSY field theories. In the last 20 years, SQM has given us deep 
insight into several aspects of standard nonrelativistic quantum mechanics. The purpose of 
these lectures is to give a brief introduction to some of these ideas. For example, 

1. It is well known that the infinite square well is one of the simplest exactly solvable 
problem in nonrelativistic QM and the energy eigenvalues are given by E n = c(n + 1) 2 
with c being a constant. Are there other potentials for which the energy eigenvalues 
have a similar form and is there a simple way of obtaining these potentials. 



2. Free particle is obviously the simplest (and in a way trivial) example in QM with no 
bound states, no reflection, and the transmission probability being unity. Are there 
some nontrivial potentials for which also there is no reflection and is it possible to 
easily construct them? 

3. Among the large number of possible potentials, only a very few are analytically solv- 
able. What is so special about these "solvable" potentials? 

4. One problem which all of us solve by two different methods (i.e. by directly solving 
Schrodinger equation and by operator method) is the one dimensional harmonic os- 
cillator potential. Can one extend this operator method to a class of potentials? In 
this context, it is worth recalling that the operator method of solving the one dimen- 
sional harmonic oscillator problem is very fundamental and in fact forms the basis of 
quantum field theory as well as many body theory. 

5. Given a potential V(x), the corresponding energy eigenvalues E n , and the scatter- 
ing matrix (i.e. the reflection and transmission coefficients R(k),T(k) in the one 
dimensional case or the phase shifts in the three dimensional case) are unique. Is 
the converse also true, i.e. given all the energy eigenvalues and R{k) and T{k) at all 
energies, is the corresponding potential unique? If not, then how does one construct 
these strictly isospectral potentials having same E n , R(k),T(k)7 

6. A related question is about the construction of the soliton solutions of the KdV and 
other nonlinear equations. Can these be easily constructed from the formalism of 
SQM? 

7. WKB is one of the celebrated semiclassical approximation scheme which is expected 
to be exact for large quantum numbers but usually is not as good for small quantum 
numbers. Is it possible to have a modified scheme which by construction would not 
only be exact for large quantum numbers but even for the ground state so that it has 
a chance to do better even when the quantum numbers are neither too large nor too 
small? Further, the lowest order WKB is exact in the case of only two potentials, 
i.e. the one dimensional oscillator and the Morse potentials. Can one construct a 
modified WKB scheme for which the lowest order approximation will be exact for a 
class of potentials? 



The purpose of these lectures is to answer some of the questions raised above. We 
shall first discuss the basic formalism of SQM and then discuss some of these issues. For 
pedagogical reason, we have kept these lectures at an elementary level. More details as 
well as discussion about several other topics can be obtained from our book [2] and Physics 
Reports |3j on this subject. A clarification is in order here. Unlike in SUSY quantum field 
theory, in SQM, super symmetric partners are not fermions and bosons, instead here, SUSY 
relates the eigenstates and S-matrix of the two partner Hamiltonians. 



2 Formalism 

Let us consider the operators 



A= _kd i = _hd 

From these two operators, we can construct two Hamiltonians E\,E 2 given by 

E x = A^A , E 2 = AA^ . (2) 

It is easily checked that 

H > = -£;i? + v ^- ViM^M-^'M, (3) 

H a = -T-7T+li(i), V 2 (x) = W\x) + -^=W'(x), (4) 
2m dx z J 2m 



The quantity W(x) is generally referred to as the superpotential in SUSY QM literature 
while Vi t2 (x) are termed as the supersymmetric partner potentials. 

As we shall see, the energy eigenvalues, the wave functions and the S-matrices of H\ 
and E 2 are related. In particular, we shall show that if E is the eigenvalue of E\ then it is 
also the eigenvalue of E 2 and vice a versa. To that end notice that the energy eigenvalues of 

(1 2) 

both E\ and E 2 are positive semi-definite (En > 0) . Further, the Schrodinger equation 
for Ei is 

Ei^ = ^(D = flW^D . ( 5 ) 
On multiplying both sides of this equation by the operator A from the left, we get 

E 2 (A^) = AA^A^) = EW(A^) . (6) 



Similarly, the Schrodinger equation for H2 

H 2 ^ = AA ^(2) = E (2)^2) (7) 

implies 

H X {AHW) = A^AAHW = E^(AH^) . (8) 

Thus we have shown that if E is an eigenvalue of the Hamiltonian H\ (H2) with eigenfunc- 
tion tp, then same E is also the eigenvalue of the Hamiltonian H2{H\) and the corresponding 
eigenfunction is Aip(A^ij)). 

The above proof breaks down in case Atp^ = 0, i.e. when the ground state is annihilated 
by the operator A (note that the spectrum of both the Hamiltonians is nonnegative) . Thus 
the exact relationship between the eigenstates of the two Hamiltonians will crucially depend 
on if Aip^ is zero or nonzero, i.e. if the ground state energy E^ is zero or nonzero. 

Aip^ ^ 0: In this case the proof goes through for all the states including the ground 
state and hence all the eigenstates of the two Hamiltonians are paired, i.e. they are related 
by (n = 0,1,2,...) 

= 4 1] > , (9) 
42) = [E CD]-i/ 2 ^(i) , (10 ) 

= im~ 1/2 A^ . (ii) 

Aijj^ = 0: In this case, E^ = and this state is unpaired while all other states of the 
two Hamiltonian are paired. It is then clear that the eigenvalues and eigenfunctions of the 
two Hamiltonians H\ and H2 are related by (n = 0, 1, 2, ...) 

E^ = 0, (12) 

1>P = [E^Y^A^l , (13) 
€1 = m~ l/2 A^ • (14) 

The equation Aip^ = can be interpreted in the following two different ways depending 
on whether the superpotential W{x) or the ground state wave function is known. In 
case W(x) is known then one can solve the equation Atp^ = by using eq. (^Q) and the 
ground state wavefunction of Hi is given, in terms of the superpotential W by 

(i) r V2m 



^=iVexp[-^— / W(y)dy}. (15) 



Instead, if i/j^ is known, then this equation gives us the superpotential W, i.e. 

W(x) = — ^ 1){X) 

Several Comments are in order at this stage. 



1. Notice that if Vn+i ( ipn 1 ) of H% (H2) is normalized then the wave function ipn 

i s a ^ so normalized. Similar remarks are obviously also valid for eigenfunctions 
in eqs. (JTUJ) and (fTTj) . 

2. The operator A (At) not only converts an eigenfunction of H\(H2) into an eigenfunc- 
tion of H2(H\) with the same energy, but it also destroys (creates) an extra node in 
the eigenfunction. 

3. Note that since the bound state wave functions must vanish at x = ±00 (or at the 
two ends) hence it is clear from eq. (|15[) that Atp^ = and A^ip^ = can never 
be satisfied simultaneously and hence only one of the two ground state energies can 
be zero. As a matter of convention we shall always choose W in such a way that 
only ift^ (if at all) is normalized so that Eq > (notice that E^ = requires that 
^(2) = 0, see eq. ©). 

4. In case since the ground state wave function of Hi is annihilated by the 
operator A, this state has no SUSY partner. Thus knowing all the eigenfunctions of 
Hi we can determine the eigenfunctions of H2 using the operator A, and vice versa 
using we can reconstruct all the eigenfunctions of Hi from those of H2 except for 
the ground state. This is illustrated in Fig. ^ 

The underlying reason for the degeneracy of the spectra of Hi and H2 can be understood 
most easily from the properties of the SUSY algebra. That is we can consider a matrix SUSY 
Hamiltonian of the form 



',(2) 



(16) 



(2) 



H 



(17) 



Hi 
. F 2J 

which contains both Hi and H2- This matrix Hamiltonian is part of a closed algebra which 
contains both bosonic and fermionic operators with commutation and anti-commutation 
relations. We consider the operators 

0' 
A 



Q 



Q 



At 




(18) 



E ( " 



Figure 1: Energy levels of two (unbroken) supersymmetric partner potentials. The action 
of the operators A and A* are displayed. The levels are degenerate except that V\ has an 
extra state at zero energy. 

in conjunction with H. The following commutation and anticommutation relations then 
describe the closed superalgebra sZ(l/l): 



[H,Q] = [H,Q*]=0, 
{Q,Q f } = H , {Q,Q} = {Q t ,Q t } = . (19) 

The fact that the supercharges Q and commute with H is responsible for the degen- 
eracy in the spectra of H\ and i?2- 

Let us now try to understand as to when is SUSY spontaneously broken and when does 
it remain unbroken. In this context, let us recall that a symmetry of the Hamiltonian (or 
Lagrangian) can be spontaneously broken if the lowest energy solution does not respect that 
symmetry, as for example in a ferromagnet, where rotational invariance of the Hamiltonian 
is broken by the ground state. We can define the ground state in our system by a two 
dimensional column vector: 



4 2) 



(20) 



Then it is easily checked that in case Aiffl / then Qi> / 0, Qty ¥= and so SUSY 
is spontaneously broken while if Atp^ = then Qtpo = 0, Q^ipo = and SUSY remains 



unbroken. Unless stated otherwise, throughout these lectures we shall be discussing the 
case when SUSY remains unbroken. 

Summarizing, we thus see that when SUSY is unbroken, then starting from an exactly 
solvable potential V(x) with n bound states and ground state energy Eq, one has V\{x) = 
V(x) — Eq whose ground state energy is therefore by construction. Using the above 
formalism, one can then immediately obtain all the n — 1 eigenstates of H%. One can now 
start from the exactly solvable Hamiltonian H2 and obtain all the n — 2 eigenstates of H%. 
In this way, by starting from an exactly solvable problem with n bound states, one can 
construct new n exactly solvable potentials H2, H3, H n+ \ with n — 1, n — 2, bound 
states. 

Illustration: Let us look at a well known potential, namely the infinite square well and 
determine its SUSY partner potential. Consider a particle of mass m in an infinite square 
well potential of width L 

V(x) = 0, < x < L , 

= 00, -00 < x < ,x > L . (21) 

The normalized ground state wave function is known to be 

4 1} = (2/L) 1 / 2 sin(vrx/L), < x < L , (22) 

and the ground state energy is 

E ~ ^ V (23) 
E °-2^- (23) 

Subtracting off the ground state energy so that the Hamiltonian can be factorized, the 

energy eigenvalues of H\ = H — Eq are 

4 1) = ^ 1 ^V,n = 0,l,2,..., (24) 

and the normalized eigenfunctions of H\ are (the same as those of H), i.e. 

V4 1} = (2/L)V2 sin (n + T 1)nX , < x < L . (25) 

The superpotential for this problem is readily obtained using eqs. (|16|) and (|2*2*)) 

W(x) = ^L-^cot(^x/L) (26) 

V2m L 




Figure 2: The infinite square well potential V = of width tt and its partner potential 
V = 2 cosec 2 x in units Ti = 2m = 1 

and hence the supersymmetric partner potential Vi is 

V *V) = Pt2\- 2 cosec 2 (7rx/L) - 1] . (27) 

The wave functions for H2 are obtained by applying the operator A to the wave functions of 
H\. In particular we find that the normalized ground and first excited state wave functions 
are 

',(2) _ m ,(2) 2 



^ = -2 J — sin 2 (vrx/L), = 7= sin(7rx/L) sin(27rx/L) . (28) 

Thus we have shown using SUSY that two rather different potentials corresponding to 
Hi and H2 have exactly the same spectra except for the fact that H2 has one fewer bound 
state. In Fig. El we show the supersymmetric partner potentials V\ and V2 and the first few 
eigenf unctions. For convenience we have chosen L = tt and fi = 2m = 1. 

We can now start from H2 and using its ground state wavefunction as given by eq. (|28|) . 
the corresponding superpotential turns out to be 

W 2 (x) = =^cot(7rx/L) , (29) 

V2m L 

and hence the partner potential Vs(x) turns out to be 

h 2 tt 2 

V 3 (x) = ——2 [6 cosecVx/L) - 4] . (30) 



Using the above superpotential W2, the ground state wavefunction of V%{x) is easily com- 
puted by using eq. (|13|) 

ip ( o\x) oc sin 3 (vr2;/L) . (31) 

In this way one can construct one (discrete) parameter family of potentials given by (p = 
0,1,2,...) 

h 2 TT 2 

V p +x(x) = ^-^[p(p + 1) cosecVx/L) - p\ (32) 
and it is easily seen that its spectrum is given by 

E^ +v > = n(n + 2p + 2) , ra = 0, 1, 2, ... , (33) 

while its eigenfunctions can be easily derived recursively from those of infinite-square well. 
For example, its ground state wave function is 

oc sm( p+1 \Trx/L). (34) 

In this way we have shown, how to generate a whole class of new solvable potentials by 
starting from an analytically solvable problem. 

So far, we have explicitly included all factors like h, m etc. However, from now onward, 
for simplicity (unless stated otherwise), we shall work in units where fi = 1, 2m = 1. 

Scattering: Supersymmetry also allows one to relate the reflection and transmission 
coefficients in situations where the two partner potentials have continuous spectra. In 
order for scattering to take place in both of the partner potentials, it is necessary that the 
potentials V\ 2 are finite as x — > —00 or as x — ► +00 or both. Let us define 

W{x -> ±00) = W± . (35) 

Then it follows that 

Vi,2 W| as nioo. (36) 

Let us consider an incident plane wave e lkx of energy E coming from the direction 
x — * —00. As a result of scattering from the potentials V\^{x) one would obtain transmitted 
waves Ti t 2(k)e lk ' x and reflected waves Ri^{k)e~ tkx . Thus we have 

V> {1,2) (A;, x -» -00) -> e ikx + R 1>2 e~ lkx , 

^(A^z-^+oo) -^T lt2 e lk ' x , (37) 



where k and k 1 are given by 

k = (E- W 2 ) 1 ' 2 , k' = {E- Wl) 1 ' 2 . (38) 

SUSY connects continuum wave functions of H\ and H2 having the same energy analogously 
to what happens in the discrete spectrum. Thus using eqs. ()13[) and ()14j) we have the 
relationships: 

e ikx + R ie - ikx = N[(-ik + W-)e ikx + (ik + W-)e- ikx R 2 ] , 

Tie ik'x = N ^_ ik ' + W+ yk'^ > (39) 

where N is an overall normalization constant. On equating terms with the same exponent 
and eliminating N, we find: 

A few remarks are in order at this stage. 

(1) Clearly \Ri\ 2 = I-R2I 2 and |Ti| 2 = IT2I 2 , that is the partner potentials have identical 
reflection and transmission probabilities. 

(2) Ri(Ti) and #2(^2) have the same poles in the complex plane except that R±(Ti) has 
an extra pole at k = —iW—. This pole is on the positive imaginary axis only if W- < in 
which case it corresponds to a zero energy bound state. 

(3) In the special case that W + = W-, we have that T±(k) = T2(k). 

(4) When W- = then R x {k) = -R 2 (k). 

(5) For symmetric potentials, W+ = —W-, and hence k = k' so that the relation between 
T\ and T2 is the same as that between R\ and R2, i.e. R\ (k) /Hq (k) = T\{k)/T2(k). 

Reflectionless Potentials: It is clear from these remarks that if one of the partner 
potentials is a constant potential (i.e. a free particle), then the other partner will be of 
necessity reflectionless. In this way we can understand the reflectionless potentials of the 
form V(x) = n(n + 1) sech 2 x which play a critical role in understanding the soliton solutions 
of the Korteweg-de Vries (KdV) equation. Let us consider the superpotential 

W(x) = B tanh x . (41) 



The two partner potentials are 

Vx = B 2 -B (B + l)sech 2 x , 

V 2 = B 2 - B(B - l)sech 2 x . (42) 

We see that for B = 1, V2(x) corresponds to a constant potential (free particle) with no 
bound states, transmission coefficient T2(k) = 1 and R2(k) = 0. Hence the corresponding 
V\ is a reflectionless potential with precisely one bound state at Eq = 0. Further by using 
eq. Q4L)j) it follows that its transmission coefficient is given by 

w-v-fiF^ry («) 

If instead we choose B = 2 then one finds that Vi = 4 — 6sech 2 x, while as seen above, 
V2 = 4 — 2sech 2 x is reflectionless potential with one bound state! Hence it follows that V± 
must be a reflectionless potential with two bound states and using the SUSY machinery (eqs. 
(|12|) to ((T3|0 . we can immediately obtain the two eigenvalues of V\ and the corresponding 
eigenf unctions. They are 

Efp (B = 2) = , ipjp oc sech 2 x , (44) 

e[ 1] (B = 2) = 3 , oc sechx tanhx. (45) 

Further using eqs. (|4()jl and ()43|) one can show that its transmission amplitude is given by 

t Or n o\ (2-ik)(l-ik) 

Tl{KB = 2) = j^W^ k)' (46) 

Generalization to arbitrary integer B is now straight forward, i.e. by choosing B = 3, 4, 
it is then easy to see that the one (discrete) parameter family of reflectionless potentials is 
given by (p = 1,2,3,...) 

Vi(x,B =p) = p(p + l)sech 2 x - p 2 , (47) 

having p bound states whose energy eigenvalues and eigenfunctions can be recursively ob- 
tained by starting from free particle and using the SUSY machinery. For example, its 
ground state energy is zero while the ground and first excited state eigenfunctions are 

iPq(x,B=p) oc sech p x , ift^(x,B=p) oc sech p ~ 1 x tanhx . (48) 



Further, its transmission coefficient is given by 

(p - ik)(p - 1 - ik)...(l - ik) T(-p - ik)T{p + 1 - ik) 
U ,P) (-p-ik)(-p + l-ik)...(-l-ik) T(-ik)T(l-ik) ' { ' 

SUSY in n-dimensions: So far we have discussed SUSY QM on the full line (—00 < 

x < 00). Many of these results have analogs for the n-dimensional potentials with spherical 

symmetry. For example, for spherically symmetric potentials in three dimensions one can 

make a partial wave expansion in terms of the wave functions: 

*Pnlm(r, 9, 0) = - r R nl {r)Y lm (0, <P) • (50) 
Then it is easily shown that the reduced radial wave function R n i satisfies (0 < r < 00) 



d 2 R n i(r) 



+ 



V(r) + 



1(1 + 1) 



R nl (r)^HR nl (r) = ER nl (r). (51) 



dr 2 

Notice that this is a Schrodinger equation for an effective one dimensional potential which 
contains the original potential plus an angular momentum barrier. 

As an illustration, let us now show how both the Coulomb and the oscillator potentials 
can be cast in SUSY formalism. 

Oscillator Potential: It is easily seen that in case we start with 

-itt / s cor l + l 

W(r) = - — , (52) 

then the two supersymmetric partner potentials are 

Vl(r) _ ^! + fci)_ ( , + 3/2) . 



4 r ^ 



Thus we see that in the oscillator case, a given partial wave I has as its SUSY partner the 
same oscillator potential but for the I + l'th partial wave. In particular, the S and P-wave 
oscillator potentials are the SUSY partners of each other. 
Coulomb Potential: Let us start with 

"v-mrr-ir- (54) 

then it is easily shown that the two supersymmetric partner potentials are 

VUr) = _^ + ll±H_l 
4(1 + 1) 2 r 2 r 

e 4 (I + !)(! + 2) e 2 



Ur) = 771 I + 



4(Z + I) 2 r 2 



(55) 



Thus we see that even for the Coulomb potential, a given partial wave I has as its SUSY 
partner the same Coulomb potential but for the I + l'th partial wave. 

As in the one dimensional case, we can also obtain the relationship between the phase- 
shifts for the (radial) partner potentials. The asymptotic form of the radial wave function 
for the l'th. partial wave is 

R(r, Q - ^[Si(k'y k ' r - {-l) l e~ lk ' r ] , (56) 

where S\ is the scattering function for the Z'th partial wave, i.e. S,(fc) = e iSl ^ and S t is the 
phase shift for the l'th partial wave. Proceeding exactly as in the one dimensional case (see 
eqs. (|35|) to (|4U|) ) . we obtain the following relationship between the scattering functions 
(and hence phase shifts) for the two partner potentials 

s ' (1, w=(raO s < 2,(fe,) ' (57) 



Here W + = W(r -> oo) and k' = - W\. 

Before finishing this section, it is amusing to note that the famous problem of charged 
particle in a uniform magnetic field can in fact be cast in the language of SQM. Actually 
one can prove a stronger result. In particular, one can show that not only uniform but even 
for an arbitrary magnetic field, the Pauli equation in two dimensions can be cast in the 
SQM formalism so long as the gyromagnetic ratio g = 2. 

Let us consider the special case when the motion of the charged particle is in a plane 
perpendicular to the magnetic field. In this case the Pauli Hamiltonian,for g = 2 has the 
form 

Hpauli = 7T- J2(Pi + - A if + 7, ^3^3 • (58) 

2m f— f c 2mc 

i=i 

One can show that in this case if we define 

Q =±{QV>+iQV>), 



6 6 

-{p y + -Ay)oi + (p x + -A x )a 2 

c c 



6 6 

(p x + -A x )ai + (p y + -A y )a 2 



(59) 

then the Hermitian supercharges Q^,Q^ and Hp au n satisfy the SUSY algebra 

{Q( a \Q®} = Hp8 ab , [H P ,Q^} = 0, a, b = 1,2. (60) 



3 Shape Invariance and Solvable Potentials 

Using the ideas of SUSY QM developed in the last section and an integrability condition 
called the shape invariance condition, we now show that the operator method for the har- 
monic oscillator can be generalized to the whole class of shape invariant potentials (SIP) 
which include essentially all the popular, analytically solvable potentials. Indeed, we shall 
see that for such potentials, the generalized operator method quickly yields all the bound 
state energy eigenvalues, eigenfunctions as well as the scattering matrix. It turns out that 
this approach is essentially equivalent to Schrodinger's method of factorization 4 although 
the language of SUSY is more appealing. 

Let us now explain precisely what one means by shape invariance. If the pair of SUSY 
partner potentials V\ )2 {x) as defined in eq. © are similar in shape and differ only in the 
parameters that appear in them, then they are said to be shape invariant. More precisely, 
if the partner potentials Vi j2 {x;a\) satisfy the condition^ 



where a\ is a set of parameters, a 2 is a function of a\ (say a 2 = f{a\)) and the remainder 
R(a±) is independent of x, then V\{x\ a\) and V 2 (x; a\) are said to be shape invariant. The 
shape invariance condition (|61[) is an integrability condition. Using this condition and the 
hierarchy of Hamiltonians discussed in the previous section, one can easily obtain the energy 
eigenvalues and the eigenfunctions of any SIP when SUSY is unbroken. 

Let us start from the SUSY partner Hamiltonians H\ and H 2 whose eigenvalues and 
eigenfunctions are related by SUSY. Further, since SUSY is unbroken we know that 



We will now show that the entire spectrum of H\ can be very easily obtained algebraically 
by using the shape invariance condition 1)61(1. As a first step, let us obtain the eigenvalue 
and the eigenfunction of the first excited state of H% . We start from eq. ()61() and on adding 
kinetic energy operator to both the sides, the shape invariance condition can also be written 
as 



Since the two Hamiltonians differ by a constant, hence it is clear that all their eigenvalues 
must differ from each other only by that constant and further all their eigenfunctions must 



V 2 (x; ai) = Vi(x; a 2 ) + -R(ai) 



(61) 




(62) 



H 2 (x; ai) = Hi(x; a 2 ) + R(ai) = A\x, a 2 )A(x, a 2 ) + R{a\) . 



(63) 



be proportional to each other. In particular 

E^im) = ^ 1) (a 2 )+i?(a 1 ), 4 2) (z,ai) « ^fo^)- (64) 

But clearly £q (a 2 ) = so long as the corresponding ground state wave function 

%'{x,a 2 ) oc exp[- / PF(t/, a 2 ) dy] , (65) 

is an acceptable ground state wave function. And hence using eqs. (j64j) and ((12)1 we obtain 
the eigenvalue of the first excited state 

E[ 1) (a 1 ) = E^(a 1 ) = R(a 1 ). (66) 

Similarly, using eqs. (|64[) and 1(13)1 we can obtain the eigenfunction of the first excited state 

ipl(x,ai) cx A\x,ai)ipQ (x,ai) oc ^(x, oi)^q ^(:c,a 2 ). (67) 

This procedure is easily generalized and we can obtain the entire bound state spectrum 
and the corresponding eigenfunctions of Hi. To that purpose, let us construct a series of 
Hamiltonians H s , s = 1, 2, 3 . In particular, following the discussion of the previous section, 
it is clear that if H\ has n bound states then one can construct n such Hamiltonians 
H2,Hs ■ ■ ■ H n+ i and the p'th Hamiltonian H p will have the same spectrum as Hi except 
that the first p — 1 levels of Hi will be absent in H p . On repeatedly using the shape 
invariance condition ((61() . it is then clear that 

H s = -^ + Vi{x-a s ) + Y^R{a k ), (68) 

where a s = / s (ai) i.e. the function / applied s — 1 times. Let us compare the spectrum 
of H s and H s+ i. In view of eqs. ((61|) and ((68)1 we have 

d 2 s 
H s+ i = + 1/1 ( x; + X/ - R ( a fc) 

= -^ + V 2 {x-a s ) + Y j R(a k ). (69) 

Thus if, and H s+ i are SUSY partner Hamiltonians and hence have identical bound state 
spectra except for the ground state of H s whose energy is 

E ( s) = J2R(a k ). (70) 
fc=l 



This follows from eq. (|68[) and the fact that Eq = 0. On going back from H s to H s _\ 
etc, we would eventually reach H 2 and H\ whose ground state energy is zero and whose 
n'th level is coincident with the ground state of the Hamiltonian H n . Hence the complete 
eigenvalue spectrum of H\ is given by 

n 

4 1) («i) = E i? (^); 4 1) («i) = °- ( 71 ) 
fc=i 

As far as the corresponding eigenfunctions are concerned, we now show that, similar 
to the case of the one dimensional harmonic oscillator, the bound state wave functions 

(x;ai) for any shape invariant potential can also be easily obtained from its ground 
state wave function ifjQ~\x; a\) which in turn is known in terms of the superpotential. This 
is possible because the operators A and link up the eigenfunctions of the same energy 
for the SUSY partner Hamiltonians Hi 2 . Let us start from the Hamiltonian H s as given 
by eq. ()68|) whose ground state eigenfunction is then given by ip^\x; a s ). On going from 
H s to H s ^i to H2 to H\ and using eq. (|14[) we then find that the n'th state unnormalized, 
energy eigenfunction ipn\x;ai) for the original Hamiltonian H\(x;a{) is given by 

ip£\x;ai) oc A\x;ai)A ] {x\a 2 )...A ] {x;a n ) ip§~\x; a n+ i) , (72) 

which is clearly a generalization of the operator method of constructing the energy eigen- 
functions for the one dimensional harmonic oscillator. 

It is often convenient to have explicit expressions for the wave functions. In that case, 
instead of using the above equation, it is far simpler to use the identity 

ip£\x;ai) = A\x;ai) ip^l 1 (x;a 2 ) . (73) 

Finally, in view of the shape invariance condition (|61|) . the relation (|4Uj) between scat- 
tering amplitudes takes a particularly simple form 

thereby relating the reflection and transmission coefficients of the same Hamiltonian Hi at 
a\ and a 2 {= f{a\)). 

Illustrative Examples: It may be worthwhile at this stage to illustrate this discussion 
with concrete examples and obtain their energy eigenvalues and eigenfunctions analytically. 



As a first illustration, let us consider the partner potentials Vx >2 (x, B) as given by eq. (|42[). 
It is easily seen that indeed these are shape invariant. In particular we observe that 

V 2 (x, B) = Vx(x, B - 1) + B 2 - (B - l) 2 , (76) 

so that in this particular case a\ = B, a 2 = B — 1 and R(a\) = B 2 — (B — l) 2 = a 2 — a 2 . 
Hence, from eq. (|71|) it immediately follows that the complete bound state spectrum of 
Vi(x, B) is given by (n = 0, 1, 2, ...) 

n 

Eg\B) = = «i - «n+i = B 2 -(B- n) 2 . (77) 

fc=i 

It is really remarkable that in two lines one is able to get the entire spectrum for this 
potential algebraically. What about the eigenfunctions? Using the superpotential as given 
by eq. (j^T|) it follows that 

^\x,B) oc sech B {x). (78) 

Clearly, this is an acceptable wave function so long as B is any positive number. On the 
other hand, the first excited state wavefunction is obtained by using the formula (|72|) 

ip^ (x, a\ = B) oc A + (x, B)i(jq(x, a 2 = B— 1) oc [— -^-+Stanhx]sech B_1 x oc sech' B_1 xtanhx 

ax 

(79) 

Proceeding in this way, one can calculate all the bound state eigenfunctions of V%(x). 

It is worth pointing out that here B is any positive number and not just an integer. 
Further it is also clear from here that the number of bound states is n + 1 in case n < B < 
n + 1. Note that only when B is integer that V± 2 are reflectionless potentials. 

As a second example, it is easily checked that the Coulomb and the oscillator potentials 
(in fact in arbitrary number of dimensions) are also examples of shape invariant potentials. 
For example, it is easy to see that in the Coulomb case, the partner potentials as given by 
eq. (|55|) are shape invariant potentials (SIP) satisfying 

V 2 (r, I, e) = Vxir, I + 1, e) + - . (80) 

4 4 

Thus in this case, a\ = I, a 2 = 1 + 1 and R(a\) = — j^p and hence the complete spectrum 
of Vi(r, I, e) is given by (n = 0, 1, 2, ...) 



Similarly, the oscillator partner potentials are SIP satisfying 



V 2 (r, l,Lo) = V 1 (r,l + l,co)-(l + l/2)u> + (Z + 5/2)w 



(82) 



so that in this case ai = l,a 2 = I + 1 while R(ai) = 2uj and hence from eq. (|71|) it 
immediately follows that the entire spectrum of the oscillator potential Vi(r, l,u) is given 



Let us now discuss the interesting but difficult question of the classification of various 
solutions to the shape invariance condition (|61|) . This is clearly an important problem 
because once such a classification is available, then one can discover new SIPs which are 
solvable by purely algebraic methods. Unfortunately, the general problem is still unsolved, 
only two classes of solutions have been found so far. In the first class, the parameters a\ 
and a 2 are related to each other by translation {a 2 = «i + ol). 

Remarkably enough, all well known analytically solvable potentials found in most text 
books on nonrelativistic quantum mechanics belong to this class. In the second class, the 
parameters a\ and a 2 are related to each other by scaling (02 = qa\). 

3.1 Solutions Involving Translation 

We shall now point out the key steps that go into the classification of SIPs in case a 2 = a\+a. 
Firstly one notices the fact that the eigenvalue spectrum of the Schrodinger equation is 
always such that the n'th eigenvalue E n for large n obeys the constraint 



where the upper bound is saturated by the infinite square well potential while the lower 
bound is saturated by the Coulomb potential. Thus, for any SIP, the structure of E n for 
large n is expected to be of the form 



by 



E^=2nto, n = 0,1,2,.... 



(83) 



A/n 2 <E n < Bn 2 



(84) 



Q 



- 2 < a < 2. 



(85) 



a 



Now, since for any SIP, E n is given by eq. I|71[). it follows that if 



R(a k )~Y, k ' 3 



(86) 



then 

- 3 < p < 1 . (87) 

How does one implement this constraint on i?(a^)? While one has no rigorous answer 
to this question, it is easily seen that a fairly general factorizable form of W(x;ai) which 
produces the above fe-dependence in R(dk) is given by 

n 

W(x; at) = J2i( k i + Ci)ft(x) + + a) + / f (s)] , (88) 

i=l 

where 

ai = (ki, k 2 ...) , a2 = (fci + a, &2 + ) 5 (89) 

with Ci,a,f3 being constants. Note that this ansatz excludes all potentials leading to E n 
which contain fractional powers of n. On using the above ansatz for W in the shape 
invariance condition eq. (|61|) one can obtain the conditions to be satisfied by the functions 
gi(x), hi(x), fi(x). One important condition is of course that only those superpotentials W 
are admissible which give a square integrable ground state wave function. 

In Tabled we give expressions for the various shape invariant potentials V\{x), superpo- 
tentials W(x), parameters a± and a 2 and the corresponding energy eigenvalues En \ Except 
for first 3 entries of this table, W(x + xq) is also a solution. 

Several remarks are in order at this time. 

1. Throughout this section we have used the convention of fi = 2m = 1. It would naively 
appear that if we had not put h = 1, then the shape invariant potentials as given in 
Tabled would all be ^-dependent. However, it is worth noting that in each and every 
case, the ^.-dependence is only in the constant multiplying the x-dependent function 
so that in each case we can always redefine the constant multiplying the function and 
obtain an ^-independent potential. For example, corresponding to the superpotential 
given by eq. (|41|). let us discuss a more general potential given by 

W(x) = Btanh(ax + x ) , (90) 

so that the corresponding h and m-dependent partner potentials are given by 

VUx, B) = W 2 - JLj^L = B 2 - B(B + -^)sech 2 (ax + x Q ) 
V2m dx V2m 

V 2 (x,A) = W 2 + — ^J?^_ = b 2 - B(B - -^L)sedi 2 (ax + x ) , (91) 



Table 1: Shape invariant potentials with (n=l,2) in which the parameters ct2 and a\ are 
related by translation (02 = a\ + (5). The energy eigenvalues and eigenfunctions are given 
in units fi = 2m = 1. The constants A, B, a, 00, 1 are all taken > 0. Unless otherwise stated, 
the range of potentials is— oo<x<oo,0<r<oo. For spherically symmetric potentials, 
the full wave function is (p n i m (r, 8, <p) = tp n l(r)Yi m (0 , <p). 



Potential 


W (X) 


Vi(x; ai) 


Ol 


Shifted oscillator 


^u)x — b 


\u 2 (x-§) 2 -^/2 


LO 


3-D oscillator 




ic 2 r 2 + ^-(/ + 3/2)c 


I 


Coulomb 


<? C+i) 


e 2 , 1(1 + 1) , e 4 


I 


2(1 + 1) r 


r -r r -z I 4(j+i)^ 


Morse 


A — B cxp (—ax) 


A 2 1 r~>2 / n \ 

A + ±( exp(— 2ax) 
-2B(^ + a/2) cxp(-ax) 


A 

A 


Scarf II 

(hyperbolic) 


Atwahax + Bsech ax 


A + (B - A 2 - Aa)sech 2 ax 
+B(2A + a)sech ax tanh ax 


A 


Rosen-Morse II 
(hyperbolic) 


.4 tanh ax + B/A 
(B < A 2 ) 


A 2 + B 2 /A 2 - A(A + a)sech 2 ax 
+ 2B tanh ax 


A 


Eckart 


-4 coth ar + B/A 
(B > A 2 ) 


A 2 + B 2 /A 2 - 2B coth ar 
+A(A — a)cosech ar 


A 


Scarf I 

(trigonometric) 


A tan ax — B sec ax 
(-iff <ax< iff) 


-A 2 + (A 2 + B 2 - Aa)sec 2 ax 
—B(2A — a) tan ax sec ax 


A 


Poschl-Teller 


^cothar - B cosech ar 
(A < B) 


A 2 + (B 2 + A 2 + Aa)cosech 2 ar 
—B(2A + a) coth ar cosech ar 


A 


Rosen-Morse I 
(trigonometric) 


-A cot ax - B/A 
(0 < ax < ff) 


A(A — a)cosec 2 ax + 2B cot ax 
-A 2 + B 2 /A 2 


A 



Note that the wave functions for the first four potentials (Hermite and Laguerre polynomials) are special 
cases of the confluent hypergeometric function while the rest (Jacobi polynomials) are special cases of the 
hypergeometric function. In the table si = s — n + a, S2 = s — n — a, S3 = a — n — s, S4 = — (s + 11 + a). 





Eigenvalue 


Variable y 




Wave function ip n (y) 


LO 


mo 


y 


= v / ¥(-- 


2b \ 

UJ j 


exp (~\y 2 )H n {y) 


l+l 


2nuo 


y 


= \^r 2 




y^ 2 M-b)Ll +1/2 (y) 


l+l 


4 4 

e e 




re 2 




/ +1 > exp(-Iy) L»+\y) 


4(i+l) 2 4(n+l + l)' A 


~ (n+i+1) 




A rv 

r\ Lx. 


A 2 — ( A — na) 2 


y 
s 


_ IB e ~ ax 
a ' 

= A/a 




y exp^ 2i/j Ll n (y) 


-T± iX 


A 2 (A -nrx\ 2 


y 

s 


= sinh ax, 
= A/a, A = 


B/a 


1 k 1 1 y ) e 

xf5 ( i A-.-l/ 2 ,-iA-.-l/2) (iy) 


A- a 


A 2 - (A - no) 2 


y 


= tanh ax, 




(l-y)^ 2 (l+y)^ 2 




~B 2 /(A~na) 2 


s 


= A/a, A = 


B/a 2 






+B 2 /A 2 


a 


= \/(s-n) 






A + a 


A 2 - (A + na) 2 


y 


= coth ar, 








-B 2 /(A + na) 2 


s 


= A/a, \ = 


B/a 2 


xPi S3 ' S4) (2/) 




+B 2 /A 2 


a 


= A/(n + s) 




A + a 


(A + na) 2 - A 2 


y 

s 


= sin ar, 
= A/a, A = 


B/a 


(!_„)(— A)/2( 1+y )(.+A)/2 
xp ( s _,_ 1/2 , s+X _ 1/2 , w 


A- a 


A 2 - (A - na) 2 


y 

s 


= cosh ar, 
= A/a, A = 


B/a 


(y-l)^>/ 2 (y+l)-^>/ 2 

xf5 (A-.-l/ a ,-A-.-l/2) (tf) 


A + a 


(A + na) 2 - A 2 


y 


= icot ax, 




(y 2 - l)-( 3 +")/ 2 exp (aax) 




-B 2 /(A + na) 2 


s 


= A/a, A = 


B/a 2 


xp (- s - n+ ia,-s-n-ia)^ 




+B 2 /A 2 


a 


= A/(s + n) 





which are shape invariant. On redefining 

n f\ 

B(B + —) = a ; (92) 
\/2m 

where a is Ti and m-independent parameter, we then have a /i-independent potential 
V{x) = V\{x) — B 2 = — asech 2 (ax + xq). 

2. It may be noted that the Coulomb as well as the harmonic oscillator potentials in 
n-dimensions are also shape invariant potentials. 

3. What we have shown here is that shape invariance is a sufficient condition for exact 
solvability. But is it also a necessary condition? The answer is clearly no. Firstly, 
it has been shown that the solvable Natanzon potentials are in general not shape 
invariant. However, for the Natanzon potentials, the energy eigenvalues and wave 
functions are known only implicitly. Secondly there are various methods one of which 
we will discuss below of finding potentials which are strictly isospectral to the SIPs. 
These are not SIPs but for all of these potentials, unlike the Natanzon case, the energy 
eigenvalues and eigenfunctions are known in a closed form. 

Solutions Involving Scaling 

From 1987 until 1993 it was believed that the only shape invariant potentials were those 
given in Tableland that there were no more shape invariant potentials. However, starting 
in 1993, a huge class of new shape invariant potentials have been discovered. It turns out 
that for many of these new shape invariant potentials, the parameters a<i and a\ are related 
by scaling {a% = qa\, < q < 1) rather than by translation, a choice motivated by the recent 
interest in g-deformed Lie algebras. Many of these potentials are reflectionless and have an 
infinite number of bound states. So far, none of these potentials have been obtained in a 
closed form but are obtained only in a series form. 

4 Strictly Isospectral Hamiltonians and SUSY 

Given any potential V(x), the corresponding bound state energy eigenvalues and the reflec- 
tion and transmission coefficients T(k), R(k) are unique. But is the converse also true? In 
particular, given the entire bound state spectrum E n and T(k),R(k) at all energies, is the 
potential uniquely defined? It turns out that the answer to this question is no. This is of 



course well known for a long time from the inverse scattering approach. In particular, it 
turns out that if a potential holds n bound states then there exist n continuous parameter 
family of strictly isospectral potentials (i.e. potentials having same energy eigenvalues and 
same reflection and transmission coefficients at all energies). These isospectral families are 
closely connected to multi-soliton solutions of nonlinear integrable systems. 

In this section we approach this problem within the SQM formalism and show that given 
any potential with at least one bound state, one can very easily construct one continuous pa- 
rameter family of strictly isospectral potentials having same E n , R(k),T(k), solongasitsgroundstatewav 
Of course this was known for a long time from inverse scattering but the Gelfand-Levitan 
approach to finding them is technically much more complicated than the supersymmetry 
approach described here. 

In SQM, the inverse scattering question can be posed as follows. Once a superpotential 
^(x) is given then the corresponding partner potentials Vi^x) are unique. But is the 
converse also true? In particular, for a given V^z) is the corresponding W(x) and hence 
Vi(x) unique? In other words, what are the various possible superpotentials W{x) other 
than W(x) satisfying 

V 2 (x) = W 2 {x) + W'(x) . (93) 

If there are new solutions, then one would obtain new potentials V\{x) = W 2 — W' which 
would be isospectral to V\ (x). To find the most general solution, let 

W(x) = W{x) + 4>(x) , (94) 

in eq. (|93j) . We then find that y(x) = (j) (x) satisfies the Bernoulli equation 

y'{x) = 1 + 2Wy , (95) 



whose solution is 



Here 



yjx) =m = ^c ln [ J ( x ) + A ] • ( 96 ) 



I{x)= [4\x')?dx\ (97) 

J — oo 

A is a constant of integration and ipQ~\x) is the normalized ground state wave function of 
V\{x) = W 2 (x) — W'(x). Thus the most general W{x) satisfying eq. (1^3^) is given by 

W(x, A) = W(x) + — ln[J(x) + A] , (98) 



so that the one parameter family of potentials 



Vi(z, A) = W 2 {x) - W'(x) = Vi{x) - In [J(x) + A] 
have the same SUSY partner V^ix). 



(99) 



The corresponding normalized ground state wave functions are 




(100) 



while the excited state wave functions are easily obtained by using relation Q14j). i.e. 



Several comments are in order at this stage. 

1. Note that this family contains the original potential V\. This corresponds to the 
choices A — ► ±oo. 

2. The fact that the potentials V\(x, A) are strictly isospectral, i.e. have same E n , R(k),T(k) 
can be seen as follows. Firstly, since all of them have the same partner potential 
V2(x), hence it follows that they must have the same spectrum. Secondly, the re- 
flection and transmission coefficients for these potentials must be related to those 
of V2(x) by the formula (|40|). But from eqs. (|HTj) and (jQHll it is easy to show that 
W(x = ±oo, A) = W± and hence the transmission and reflection coefficients for the 
entire one parameter family of potentials are identical. 

3. Note that while E n , R(k),T(k) are the same for the entire family of potentials, the 
wave functions for all these potentials are different. As is well known from inverse 
scattering, a potential is uniquely fixed only when one specifies E n , R(k),T(k) and 
normalization constant of one of the eigenfunction. 

4. Since I varies from to 1 hence if we choose the continuous parameter A so that either 
A > or A < - 1 then we ensure that the corresponding ground state wavefunctions 
tpi(x,X) are nonsingular. 

5. The one parameter family of strictly isospectral potentials Vi(x, A) can be shown to 
satisfy infinite number of conserved charges Qi (i=l,2,3,...). Two of these are 




(101) 




(102) 



which satisfy %^ = as can be easily checked by using eqs. (j98|) to The 
conserved charge Q2 implies that the one continuous parameter family of potentials 
are such that if one plots these potentials as a function of x, then the area under 
the curve is the same for the entire family of potentials. Similarly, all the other 
conservation laws tighten these potentials. 

To elucidate this discussion, it may be worthwhile to explicitly construct the one- 
parameter family of strictly isospectral potentials corresponding to the one dimensional 
harmonic oscillator. In this case 

W(x) = ~x , (103) 



so that 



LO 2 2 



Vi(x) = —x l - - . (104) 



The normalized ground state eigenfunction of V\{x) is 

^\x) = (|-) 1/4 exp(-^r 2 /4) . (105) 

Using eq. (|97|) it is now straightforward to compute the corresponding Z{x). We get 

Z{x) = 1 - - erfc l—x) ; erfc(x) = J_ T e ~t 2 dt _ ( 106 ) 

2 y 2 j v 71 " Jx 

Using eqs. (|99|) and ()100|) . one obtains the one parameter family of isospectral potentials 
and the corresponding wave functions. In Figs. |3] and ^ , we have plotted some of the 
potentials and the corresponding ground state wave functions for the case to = 2. 

We see that as A decreases from 00 to 0, V\ starts developing a minimum which shifts 
towards x = —00. Note that as Ai finally becomes zero this attractive potential well is lost 
and we lose a bound state. The remaining potential is called the Pursey potential Vp(x). 
The general formula for Vp(x) is obtained by putting A = in eq. (|99|) . An analogous 
situation occurs in the limit A = — 1, the remaining potential being the Abraham-Moses 
potential. 



5 Supersymmetry Inspired WKB Approximation 

We will describe here a recent extension of the semiclassical approach inspired by super- 
symmetry called the super symmetric WKB (SWKB) method. It turns out that for many 
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Figure 3: Selected members of the family of potentials with energy spectra identical to the 
one dimensional harmonic oscillator with u = 2. The choice of units is H = 2m = 1. The 
curves are labeled by the value of A, and cover the range < A < oo. The curve A = oo is 
the one dimensional harmonic oscillator. The curve marked A = is known as the Pursey 
potential and has one bound state less than the oscillator. 




Figure 4: Ground state wave functions for all the potentials shown in Fig. EJ except the 
Pursey potential. 



problems the SWKB method gives better accuracy than the WKB method. Further, we 
discuss and prove the remarkable result that the lowest order SWKB approximation gives 
exact energy eigenvalues for all SIPs with translation which include essentially all analyt- 
ically solvable potentials discussed in most text books on QM. In this section we do not 
assume h = 2m = 1 but explicitly keep all the factors of h and m. 

The semiclassical WKB approximation for one dimensional potentials with two classical 
turning points is discussed in most quantum mechanics textbooks. The lowest order WKB 
quantization condition is given by (n = 0, 1,2, ...) 

f R dxJ2m[E - V(x)} = (n + l/2)ftvr . (107) 

J XL 

In the special case of the one dimensional harmonic oscillator and the Morse potential, it 
turns out that the lowest order WKB approximation ()107|) is in fact exact and further, the 
higher order corrections are all zero. 

About twenty years ago, combining the ideas of SUSY with the lowest order WKB 
method, the lowest order SWKB quantization condition was obtained in case SUSY is 
unbroken and it was shown that it yields energy eigenvalues which are not only accurate for 
large quantum numbers n but which, by construction, are also exact for the ground state 
(n = 0). We shall now discuss the SWKB quantization in detail. 

For the potential V\(x) corresponding to the superpotential W(x), the lowest order 
WKB quantization condition (|107j) takes the form 



■>R 



1 2m 



EjP - W 2 (x) + -JL=W'(x) 
V2m 



dx = (n + l/2)hir. (108) 



Let us assume that the superpotential W(x) is formally 0(h ). Then, the W term is clearly 
0(h). Therefore, expanding the left hand side in powers of ft gives 

f \/2m[E^ -W 2 (x)\ dx + ^- [ dx = + ... = (n + 1/2)^, (109) 

Ja 2Ja ^E { n ] -W 2 {x) 

where a and b are the turning points defined by En^ = W 2 (a) = W 2 (b). The 0(h) term in 
eq. (|109|) can be integrated easily to yield 



& • -i 
— sin 

2 



W(x) 



(110) 



In the case of unbroken SUSY the superpotential W(x) has opposite signs at the two turning 
points, that is 



-W(a) = W(b) = \/E ( n ) . (Ill) 

For this case, the 0(h) term in ()110|) exactly gives Htt/2, so that to leading order in fi the 
SWKB quantization condition when SUSY is unbroken is 

f ^2m[E ( n 1] - W 2 {x)} dx = nhir, n = 0,1,2,.... (112) 

J a 

Proceeding in the same way, the SWKB quantization condition for the potential V^x) turns 
out to be 

f \j2m[E^ - W 2 (x)} dx = (n + n = 0,1,2,.... (113) 

J a 

Some remarks are in order at this stage. 

(i) For n = 0, the turning points a and b in eq. (|112|) are coincident and = 0. Hence 
the SWKB condition is exact by construction for the ground state energy of the potential 
V x {x). 

(ii) On comparing eqs. (|112j) and (j!13|) . it follows that the lowest order SWKB quantiza- 
tion condition preserves the SUSY level degeneracy i.e. the approximate energy eigenvalues 
computed from the SWKB quantization conditions for V\(x) and V^a;) satisfy the exact 
degeneracy relation E^ x = E$ ■ 

(iii) Since the lowest order SWKB approximation is not only exact, as expected, for large 
n, but is also exact by construction for n = 0, hence, unlike the ordinary WKB approach, 
the SWKB eigenvalues are constrained to be accurate at both ends, at least when the 
spectrum is purely discrete. One can thus reasonably expect better results than the WKB 
scheme even when n is neither small nor very large. 

(iv) For spherically symmetric potentials, unlike the conventional WKB approach, in the 
SWKB case one obtains the correct threshold behaviour without making any Langer-like 
correction. This happens because, in this approach 

5*0 ~ [E-W 2 ) 1 ' 2 r ~° -%h(l + l)/r, (114) 

so that 



ip(r) ~ exp 



r 



Sodr 



r ~° r l+1 . (115) 



5.1 Exactness of the SWKB Condition for Shape Invariant Potentials 

In order to determine the accuracy of the SWKB quantization condition as given by eq. 
(|112|) . researchers first obtained the SWKB bound state spectra of several analytically 
solvable potentials. Remarkably they found that the lowest order SWKB condition gives 
the exact eigenvalues for all SIPs with translation! Let us now prove this result. 
Recall that the shape invariance condition eq. (|6Tj) on the partner potentials is 

V 2 {x,ai) = Vi(x,a 2 ) + R(ai) , 

where a\ is a set of parameters, a 2 is a function of a\ (say a 2 = /(ai)) and the remainder 
R(ai) is independent of x. 

In Section 3, we showed using factorization and the Hamiltonian hierarchy that the 
general expression for the s'th Hamiltonian was given by 

s-1 

4- 1/; ( -r n.l-l. 

2m dx 2 



fr 2 d 2 



fc=i 

where a s = f s ^ 1 (ai) i.e. the function / applied s — 1 times. 

The proof of the exactness of the bound state spectrum eq. (|71[) i n the lowest order 
SWKB approximation now follows from the fact that the SWKB condition (|112j) preserves 
(a) the level degeneracy and (b) a vanishing ground state energy eigenvalue. For the hier- 
archy of Hamiltonians as given by eq. ()68|) . the SWKB quantization condition takes 
the form 

r r fz. 1 1 

dx = nhir . (116) 
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Now, since the SWKB quantization condition is exact for the ground state energy when 
SUSY is unbroken, hence 

4 s) = |>(a fe ) (117) 

k=l 

must be exact for Hamiltonian as given by eq. H116|) . One can now go back in sequential 
manner from H(") to H^-V to a nd and use the fact that the SWKB method 
preserves the level degeneracy E^+i = En . On using this relation n times, we find that 
for all SIPs, the lowest order SWKB condition gives the exact energy eigenvalues. 

This is a substantial improvement over the usual WKB formula eq. ()107|) which is not 
exact for most SIPs. Of course, one can artificially restore exactness by ad hoc Langer-like 



corrections. However, such modifications are unmotivated and have different forms for dif- 
ferent potentials. Besides, even with such corrections, the higher order WKB contributions 
are non-zero for most of these potentials. 

What about the higher order SWKB contributions? Since the lowest order SWKB 
energies are exact for shape invariant potentials, it would be nice to check that higher order 
corrections vanish order by order in H. By starting from the higher order WKB formalism, 
one can readily develop the higher order SWKB formalism. It has been explicitly checked 
for all known SIPs (with translation) that up to 0(H 6 ) there are indeed no corrections.. 
This result can be extended to all orders in fi. 

Let us now compare the merits of the WKB and SWKB methods. For potentials for 
which the ground state wave function (and hence the superpotential W) is not known, 
clearly the WKB approach is preferable, since one cannot directly make use of the SWKB 
quantization condition Q112|) . On the other hand, we have already seen that for shape in- 
variant potentials, SWKB is clearly superior. An obvious interesting question is to compare 
WKB and SWKB for potentials which are not shape invariant but for whom the ground 
state wave function is known. An extensive study of several potentials indicate that by and 
large, SWKB does better than WKB in case the ground state wave function and hence the 
superpotential W is known. These studies also support the conjecture that shape invariance 
is perhaps a necessary condition so that the lowest order SWKB reproduce the exact bound 
state spectrum. 

6 New Periodic Potentials Prom Super symmetry 

So far we have considered potentials which have discrete or discrete plus continuous spectra 
and by using SUSY QM methods we have generated new solvable potentials. In this section 
we extend this discussion to periodic potentials and their band spectra. The importance of 
this problem can hardly be overemphasized. For example, the energy spectrum of electrons 
on a lattice is of central importance in condensed matter physics. In particular, knowledge 
of the existence and locations of band edges and band gaps determines many physical 
properties of these systems. Unfortunately, even in one dimension, there are very few 
analytically solvable periodic potential problems. We show in this section that SQM allows 
us to enlarge this class of solvable periodic potential problems. 



We start from the Hamiltonians H\ 2 in which the SUSY partner potentials V\ 2 are 
periodic nonsingular potentials with a period L. In view of the periodicity, one seeks 
solutions of the Schrodinger equation subject to the Bloch condition 

ip(x + L) = e ikL ?p(x) , (118) 

where k is real and denotes the crystal momentum. As a result, the spectrum shows energy 
bands whose edges correspond to kL = 0, ir, that is the wave function at the band edges 
satisfy tp{x+L) = ±ip(x). For periodic potentials, the band edge energies and wave functions 
are often called eigenvalues and eigenfunctions, and we will use this terminology in these 
lectures. In particular the ground state eigenvalue and eigenfunction refers to the bottom 
edge of the lowest energy band. 

Let us first discuss the question of SUSY breaking for periodic potentials. Since Hi = 
A^A and H% = AA^ are formally positive operators their spectrum is nonnegative and 
almost the same. The caveat "almost" is needed because the mapping between the positive 
energy states of the two does not apply to zero energy states. 

The Schrodinger equation for H\ t 2 has zero energy modes given by 

4 1,2) (x) = exp f =f J X dyW(y)) , (119) 
provided ip^ ' belong to the Hilbert space. Supersymmetry is unbroken if at least one of 

(1 2) 

the ipQ is a true zero mode while otherwise it is dynamically broken. 

(1 2) 

For a non-periodic potential we have seen that at most one of the functions ipQ ' can 
be normalizable and hence an acceptable eigenfunction. By convention we are choosing W 
such that only Hi (if at all) has a zero mode. 

Let us now consider the case when W (and hence V12) are periodic with period L. 
Now the eigenfunctions including the ground state wave function must satisfy the Bloch 
condition (|118[) . But, in view of eq. (|119|) we have 

4 1,2) (x + L) = e ±<t>L Vo 1,2) ( x ) > ( 12 °) 

where 

<t> L = f X+ W(y)dy. (121) 

J X 

(I 2) 

On comparing eqs. (|118|) and (|12U|) it is clear that for either of the wave functions V'o ' 
to belong to the Hilbert space, we must identify ±0l = ikL. But (f>L is real (since W and 



W(y)dy = (122) 



hence V\p are assumed to be real), which means that </>£, = 0. Thus, the two functions 

(1 2) 

' either both belong to the Hilbert space, in which case they are strictly periodic with 

(1 2) (1 2) 

period L: tpQ ' '{x + L) = ip ' (x), or (when <j>L / 0) neither of them belongs to the Hilbert 
space. Thus in the periodic case, irrespective of whether SUSY is broken or unbroken, the 
spectra of V\$ is always strictly isospectral. 
To summarize, we see that 

4>l= [ L 

JO 

is a necessary condition for unbroken SUSY, and when this condition is satisfied then H\p. 
have identical spectra, including zero modes. In this case, using the known eigenfunctions 
^n\x) Of Vl (x) one can immediately write down the corresponding (un-normalized) eigen- 

(2) 

functions ipn (x) of V^x). In particular, from eq. (7.3) the ground state of V^(x) is given 
by 

4% )= * = e r^ , (123) 

1>o \ x ) 

while the excited states tpn^ (x) are obtained from ipn^ (x) by using the relation 

^ ){X) = [ ~i + W ^n\^) ,(">!)• (124) 

Thus by starting from an exactly solvable periodic potential Vi(x), one gets another strictly 
isospectral periodic potential V^x). 

At this stage, it is worth pointing out that there are some special classes of periodic 
superpotentials which trivially satisfy the condition (J122j) and hence for them SUSY is 
unbroken. For example, suppose the superpotential is antisymmetric on a half-period: 

W(x + ^) = -W(x). (125) 



Then, 



Vx, 2 (x + |) = W 2 (x + ^) T W(x + ~) = U 2 ,i(x) . (126) 



Thus in this case Vi^ are simply translations of one another by half a period, and hence 
are essentially identical in shape. Therefore, they must support exactly the same spectrum, 
as SUSY indeed tells us they do. Such a pair of isospectral Vi,2 that are identical in 
shape are termed as "self-isospectral" . A simple example of a superpotential of this type is 
W(x) = cosx, so that V2(x) = cos 2 x — sinx = Ui(x + 7r). In a way, self-isospectral potentials 
are uninteresting since in this case, SUSY will give us nothing new. 



More generally, if a pair of periodic partner potentials V\2 are such that V 2 is just the 
partner potential V% up to a discrete transformation-a translation by any constant amount, 
a reflection, or both, then such a pair of partner potentials are termed as "self-isospectral 
". For example, consider periodic superpotentials that are even functions of x: 

W{-x) = W(x), (127) 

but which also satisfy the condition (|122j) . Since the function dW{x)/dx is now odd hence 
it follows that 

Vi, 2 (-x) = V 2 , 1 (x). (128) 

The partner potentials are then simply reflections of one another. They therefore have 
the same shape and hence give rise to exactly the same spectrum. A simple example of a 
superpotential of this type is again W(x) = cosx, so that V 2 {x) = cos 2 :r — since = V\{—x). 

It must be made clear here that not all periodic partner potentials are self-isospectral 
even though they are strictly isospectral. Consider for example, periodic superpotentials 
that are odd functions of x: 

W(-x) = -W{x). (129) 

Then the condition ()122j) is satisfied trivially and hence SUSY is unbroken. The function 
dW(x)/dx is even and thus ^1,2(2;) are also even. In this case, V±(x) are not necessarily 
related by simple translations or reflections. For example, the superpotential W^cc) = 
A since + B sin2cc gives rise to an isospectral pair which is not self-isospectral. On the other 
hand, W{x) = A since + B sin3cc gives rise to a self-isospectral pair since this W satisfies 
the condition (J125|) . 

6.1 Lame Potentials and Their Supersymmetric Partners 

The classic text book example of a periodic potential which is often used to demonstrate 

band structure is the Kronig-Penney model , 

00 

V(x) = Vq5{x - nL) . (130) 

—00 

It should be noted that the band edges of this model can only be computed by solving a 
transcendental equation. 

Another well studied class of periodic problems consists of the Lame potentials 

V(x,m) = p m sn 2 (x,m) , p = a(a + 1) . (131) 



Here sn(x,m) is a Jacobi elliptic function of real elliptic modulus parameter m(0 < m < 1) 
with period 4K(m), where K(m) is the " real elliptic quarter period " given by 

K(m) = / f d ° . . (132) 

Jo v 1-msir S 

For simplicity, from now on, we will not explicitly display the modulus parameter m as 
an argument of Jacobi elliptic functions unless necessary. Note that the elliptic function 
potentials (|131|) have period 2K(m). They will be referred to as Lame potentials, since 
the corresponding Schrodinger equation is called the Lame equation in the mathematics 
literature. It is known that for any integer value a = 1,2,3,..., the corresponding Lame 
potential has a bound bands followed by a continuum band. All the band edge energies 
and the corresponding wave functions are analytically known. We shall now apply the 
formalism of SUSY QM and calculate the SUSY partner potentials corresponding to the 
Lame potentials as given by eq. Q131|) and show that even though a = 1 Lame partners are 
self-isospectral, for a > 2 they are not self-isospectral. Consequently, SUSY QM generates 
new exactly solvable periodic problems! 

Before we start our discussion, it is worth mentioning a few basic properties of the Jacobi 
elliptic functions snx, cnx and dnx which we shall be using in this discussion. First of all, 
whereas sni and cnx have period 4K(m), dnx has period 2K(m) [i.e. dn(x + 2K{m)) = 
dnx]. They are related to each other by 

m sn 2 x = m — m cn 2 x = 1 — dn 2 x . (133) 

Further, 

d d d 

— snx = cnx dnx ; —cnx = — snx dnx , —dnx = — m snx cnx . (134) 
dx dx dx 

Besides 



, cnx . ts~\ n snx i / \/\ — m ,,„_<. 

sn(x + K) = ; cn(x + K) = -VI - m- ; dn(x + K) = . (135) 

dn x dn x dn x 

Finally, for m = 1(0), these functions reduce to the familiar hyperbolic (trigonometric) 
functions, i.e. 



sn(x, m = 1) = tanh x ; cn(x, m = 1) = sech x ; dn(x, m = 1) = sech x , 

sn(x, m = 0) = sin x ; cn(x, m = 0) = cos x ; dn(x, m = 0) = 1 . (136) 



It may be noted that when m = 1, the Lame potentials ()131|) reduce to the well known 
Poschl-Teller potentials 

V(x, m = 1) = a(a + 1) — a(a + l)sech 2 x , (137) 

which for integer a are reflectionless and have a bound states. It is worth adding here that in 
the limit m — > 1, K(m) tends to oo and the periodic nature of the potential is obscure. On 
the other hand, when m = 0, the Lame potential ()131[) vanishes and one has a rigid rotator 
problem (of period 2K(m = 0) = tt), whose energy eigenvalues are at E = 0, 1,4,9, ... with 
all the nonzero energy eigenvalues being two-fold degenerate. 

Finally, it may be noted that the Schrodinger equation for finding the eigenstates for 
an arbitrary periodic potential is called Hill's equation in the mathematical literature. A 
general property of the Hill's equation is the oscillation theorem which states that for a 
potential with period L, the band edge wave functions arranged in order of increasing 
energy E < E 1 < E 2 < E 3 < E A < E 5 < E e < ... are of period L, 2L, 2L, L, L, 2L, 2L, ... . 
The corresponding number of (wave function) nodes in the interval L are 0, 1, 1, 2, 2, 3, 3, ... 
and the energy band gaps are given by Ai = E% — E\, A2 = £4 — -E3, A3 = Eq — E§, ... . 
We shall see that the expected m = limit and the oscillation theorem are very useful in 
making sure that all band edge eigenstates have been properly determined or if some have 
been missed. 

Let us first consider the Lame potential ()131|) with a = 1 and show that in this case 
the SUSY partner potentials are self-isospectral. The Schrodinger equation for the Lame 
potential with a = 1 can be solved exactly and it is well known that in this case the spectrum 
consists of a single bound band followed by a continuum band. In particular, the eigenstates 
for the lower and upper edge of the bound band are given by 

Eo = m; ^o(^) = dnx, (138) 

E 1 = l; ip 1 (x) = cnx. (139) 
On the other hand, the eigenstate for the lower edge of the continuum band is given by 

E 2 = l + m; ^j 2 (x)=snx. (140) 

i.e. it extends from 1 + m to 00. Note that at m = the energy eigenvalues are at 0, 1 as 
expected for a rigid rotator and as m — > 1, one gets V(x) — > 2 — 2 sech 2 x, the band width 
1 — m vanishes as expected, and one has an energy level at E = 1. 



Using eq. (|138j) the corresponding superpotential turns out to be 

jt t / \ snx cm 

W(x) = m — : , (141) 

dnx 

On making use of eq. (|135|) it is easily shown that this W satisfies the condition (|125|) and 
hence the corresponding partner potentials are indeed self-isospectral. 

The Lame potential (j!31j) with a = 1, is one of the rare periodic potentials for which 
the dispersion relation between E and crystal momentum k is known in a closed form. 

In view of this result for a = 1, one might think that even for higher integer values of a, 
the two partner potentials would be self-isospectral. However, this is not so, and in fact for 
any integer a(> 2), we obtain new exactly solvable periodic potential. As an illustration, 
consider the Lame potential (j!31j) with a = 2. For the a = 2 case, the Lame potential has 
2 bound bands and a continuum band. The energies and wave functions of the five band 
edges are well known. The lowest energy band ranges from 2 + 2m — 25 to 1 + m, the second 
energy band ranges from 1 + 4m to 4 + m and the continuum starts at energy 2 + 2m + 25, 
where 5 = \/l — m + m 2 . 

Note that in the interval 2K{m) corresponding to the period of the Lame potential, the 
number of nodes increase with energy. In order to use the SUSY QM formalism, we must 
shift the Lame potential by a constant to ensure that the ground state (i.e. the lower edge 
of the lowest band) has energy E = 0. As a result, the potential 

V 1 (x) = -2-2m + 25 + 6msn 2 x, (142) 

has its ground state energy at zero with the corresponding un-normalized wave function 

(x) = I + m + 5 — 3m sn 2 x . (143) 

The corresponding superpotential is 

d 6msnx cnxdnx 
W = - Tx lo g ^(x) = , (111, 



and hence the partner potential corresponding to (|142[) is 



72m 2 sn 2 x cn 2 x dn 2 x 
[1 + m + 5 — 3m sn 2 x] 



U 2 (x) = -U 1 (x)+ ri , _ , , o _ . (145) 



Although the SUSY QM formalism guarantees that the potentials V\ 2 are isospectral, 
they are not self-isospectral, since they do not satisfy eq. (j!26|) . Therefore, V^ix) as given 



by eq. (|145j) is a new periodic potential which is strictly isospectral to the potential (|142j) 
and hence it also has 2 bound bands and a continuum band. Similar conclusions are also 
valid for Lame potentials (|131f) with a = 3,4,5, .... 
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